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Abstract. 8000 images of the Solar corona were captured during the June 2001 
total Solar eclipse. New software for the alignment of the images and an automated 
technique for detecting intensity oscillations using multi scale wavelet analysis were 
developed. Large areas of the images covered by the Moon and the upper corona 
were scanned for oscillations and the statistical properties of the atmospheric effects 
^ ■ were determined. The a Trous wavelet transform was used for noise reduction and 

Monte Carlo analysis as a significance test of the detections. The effectiveness of 
00 . those techniques is discussed in detail. 
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1. Introduction 



5_i ■ A number of theoretical predictions exist for the propagation of magne- 

"£o . tohydrodynamic (MHD) waves associated with coronal loops. Roberts 

et al. (1984) studied the various types of propagation through a low-/? 
plasma using reasonable approximations for the conditions inside coro- 
^ | nal loops. Their work was followed by a large number of authors since 

5_i ' then (for example see one of the many review papers by Nakariakov, 

2003) and was recently confirmed and refined by applying numerical 
modelling Nakariakov et al. (2004). 

Many attempts to observe propagations in the coronal loops have 
been made since the first theoretical predictions. One of the most chal- 
lenging types of oscillations to observe are the fast sausage mode MHD 
that have expected periodicities below lmin (see Aschwanden, 2004 for 
a detailed review). Radio, optical and X-rays observations have been 
used to detect such waves with limited success. In this paper we will 
present the application of image processing techniques as a way to 
enhance optical observations made by the Solar Eclipse Coronal Imag- 
ing System (SECIS) project during the June 2001 total solar eclipse 
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in Zambia. A detailed description of the instrument can be found in 
Phillips et al. (2000). 

Starting with Koutchmy et al. (1983) a number of authors have 
published possible detections of oscillations with periodicities below 
lOsec. Pasachoff & Landman, Pasachoff & Ladd (1984, 1987) have re- 
ported possible detections of optical intensity oscillations with periods 
in the range of 0.5-4sec, while more recently Williams et al., Williams 
et al, Katsiyannis et al. (2001, 2002, 2003) (hereafter W01, W02 and 
K03 respectively) provided strong indications of oscillations with pe- 
riodicities ~6sec while reporting on optical SECIS August 1999 total 
solar eclipse observations. 

Continuing the work published for the SECIS 1999 observations, 
we have analysed observations made during the June 2001 total solar 
eclipse in Lusaka, Zambia. Based on experience from the analysis of 
the 1999 data set, a number of numerical techniques were used in order 
to improve the signal-to-noise (S/N) ratio and establish an objective, 
numerical criterion for the identification of the corona intensity oscilla- 
tions over any statistical effects. A brief description of the observations 
and data analysis is presented here with more emphasis given to the 
advanced mathematic techniques used in an effort to improve the S/N 
ratio and determine the "real" detections of corona loop waves over the 
influences of noise in the data set. 



2. Observations & Data Reduction 

A detailed description of the SECIS instrument, as used for the 1999 
observations is provided by Phillips et al. (2000), while a discussion 
of the improvements made for the 2001 observations can be found in 
Katsiyannis et al. (2004) (hereafter K04). The observations taken by 
SECIS in 2001 and their data reduction will not be presented in detail 
in this paper as they are the subject of K04. However, a brief description 
of the data set is to follow as needed for the presentation of the image 
processing techniques reported here. 

8000 Fe xiv images of 512x512 pixel 2 with a resolution of ps 4 arc- 
sec pixel -1 were taken during the ~ 3.5 min of totality. Although the 
observing field was large enough to include the whole Moon disk and the 
lower part of the corona, we chose only to observe the North-East limb. 
This decision is in line with the 1999 observations and was taken to 
avoid edge effects of the CCD and optics as well as to include important 
parts of the outer corona. 

A brief description of the data reduction of the eclipse 2001 observa- 
tions is included here for the purpose of describing the data that were 
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used for the application of the a Trous wavelet transform and Monte 
Carlo analysis. The images taken during the 2001 observations were 
reduced by using dark and flat-field frames , for current subtraction 
and flat field correction. The 8000 images were then automatically co- 
aligned using the edge of the Moon as a reference point for a first order 
alignment. For this first alignment the moon was effectively considered 
stationary during the duration of the eclipse. A more accurate align- 
ment was subsequently achieved by using a clear feature from an area 
of the lower corona as reference. This second alignment corrected for 
the motion of the moon in respect to the solar corona during totality. 
K04 provides a detailed discussion on the alignment technique used and 
its various steps. 

After co-alignment the 8000 images of the observations form a three 
dimensional data array. The basic technique used for the detection of 
intensity oscillations throughout the SECIS project is the continuous 
wavelet transformation of the time series that corresponds to each of 
pixels of the aligned images. Details on the transformation function and 
its implementation in time series can be found in Torrence & Compo 
(1998), while examples of the application of this analysis to SECIS data 
can be found in a number of publications (e.g. W01, W02, K03 and 
K04). Additionally K03 explicitly mentions a number of criteria used 
for a wavelet detection to be considered as a solar intensity oscillation 
(as oppose to detections created by noise). To test the performance of 
K03's criteria, K04 used automated software to scanned large areas of 
the image covered by the Moon and upper corona and confirmed their 
effectiveness. 



3. Noise Reduction using the it A Trous wavelet algorithm 

3.1. S/N RATIO LIMITATIONS ON SECIS ECLIPSE DATA 

On of the most significant limitations of the SECIS eclipse observations 
is the low S/N ratio. Although the observations were taken using a 
broad Fe xiv filter and the solar corona is known to be bright in Fe xiv 
emission, there are three major factors that severely limit the S/N 
ratio achieved by SECIS. For the purposes of this paper we will only 
emphasise the S/N ratio limiting factors: 

1. The prime mirror of the SECIS telescope has a 200 mm diameter 
and a focal length of f/10. This is because the instrument was 
designed mainly with observing solar eclipses in mind and has to 
be lightweight and easy to travel. 
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2. The CCD cameras took observations at a ratio of ^40 frames per 
second. This has the obvious disadvantages of a very short exposure 
time and a very fast readout speed. With current technology such 
fast CCD readout speeds increase the readout noise drastically. 

3. As the purpose of the SECIS project is to detect high frequency 
intensity oscillations, the atmospheric effects became significant. 
This is because the earth's atmosphere itself is known to oscillate 
in these frequencies and that is causing non-Gaussian noise on the 
data set. 

3.2. Advantages of using the a Trous wavelet transform in 

OUR DATA SET. 

Having the above limitations in mind, the a Trous filtering was in- 
vestigated as a means of noise reduction because the algorithm has a 
number of advantages Starck & Murtagh (2002). 

— The computational requirements are within acceptable levels. This 
is particularly important as both of the SECIS eclipse observations 
will have a size of more than 8 GB. 

— The reconstruction algorithm is trivial. This is important as it 
makes the reconstruction of the time series more accurate. 

— The transform is known for every sample of the time series of every 
pixel. This is important to this project as the exact moment an 
oscillation starts and ends on a given pixel can be very significant. 
W02 made some interesting measurement of the propagation speed 
of a travelling wave on SECIS 1999 total solar eclipse observations 
by determining the exact time the oscillation arrived at any given 
pixel. 

— The transform is clearly evolving through the different scale in a 
predictable mater. This makes easier to choose the right scale for 
the filtering of a certain data set. 

— The transform is isotropic. As the SECIS data are also isotropic 
(specially in the time domain), any filtering used should also be 
isotropic to avoid an artifacts being introduced. 

The a Trous filtering is a relatively recent (some of the first appli- 
cation were described by Holschneider et al., 1989), s sophisticated, 
highly tunable and complicated multiresolution algorithm. Due to its 
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Figure 1. The 2001 total solar eclipse as observed by SECIS. This image is produced 
by averaging over the time axis the three-dimensional aligned array. After alignment 
the Sun remains in the same pixel area and the Moon crescent moves in respect 
to the solar corona. Areas of the Moon, outer corona and lower corona close to 
NOAA Active Region 9513 (AR 9513) are highlighted. These three were the areas 
used to test the effectiveness of the a Trous noise filtering and the Monte Carlo 
randomisation analysis and were also the exact areas also used by K04. 

complication the exact description of the algorithm is outside the pur- 
poses of this paper. More information on the algorithm, its various 
parameters, advantages and disadvantages of the various procedures 
that can be used in conjunction with the a Trous wavelet transform 
and examples of its application on astrophysical data can be found in 
Starck & Murtagh (2002). 

3.3. The application of a Trous noise reduction algorithm 

TO SECIS DATA. 

Each pixel of the three-dimensional data array of the reduced SECIS 
2001 observations was treated as a independent time series and was 
transformed to a number of coefficients on a multi-scale domain using 
the a Trous wavelet transform. B3 splines were used for correlation 
and the noise was assumed to be Gaussian. The sigma of the noise 
in the coefficients was determined automatically and multiresolution 
hard K-Sigma thresholding was used to remove the coefficients that 
were found to be noise. The time series was then reconstructed using 
the coefficients corrected for noise. 
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Figure 2. A 30x30 pixels Moon area of the 2001 observations. All pixels in that area 
had their Gaussian noise removed by applying the d Trous wavelet transformation. 
Marked with 'x' are the 374 detections made by the automated software of K04 after 
the removal of Gaussian noise. 



To test the effectiveness of the a Trous wavelet transformation in 
noise reduction the same three parts of the data set that were chosen 
by K04 were used again. Figure 1 contains the time-average image of 
the aligned SECIS 2001 total eclipse observations. After alignment the 
three-dimensional data array was averaged over the time axis resulting 
to a two-dimensional image. The edges of the image (one hundred pixels 
of the left, right and bottom of the image and two hundred pixels 
from the top) were discarded as they suffer from edge effects of the 
CCD. Highlighted are three areas of the data set that were used to 
test the effectiveness of the a Trous filtering algorithm and the Monte 
Carlo randomisation test. These are areas of the image covered by 
the Moon's disk, the outer corona and lower corona. Those areas were 
chosen to be in the proximity of AR 9513 as this is potentially the 
most interesting (from the coronal loop oscillations point of view) area 
of the 2001 observations. The proximity of the test areas to the area 
were more detections of coronal waves are expected is important as it 
produces more reliable statistics as there are no effects from large-scale 
variation in CCD sensitivity or deferent atmospheric conditions. 

The same automated detection algorithm was used as in K04 and 
the results for the periodicity range of 7-8 sec are displayed in Figures 
2, 3 and 4. Figure 2 contains 374 pixels that have false (i.e. not caused 
by the solar corona) detections of oscillations out of the 900 pixels of 
the sample. Before the d Trous noise filtering the same area contained 
5 oscillating pixels. Figure 3 contains 1054 oscillating pixels out of 
2500 while before the application of the d Trous wavelet transformation 
algorithm there were 11. On figure 4 we have detections on 276 pixels 
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Figure 3. A 50x50 pixels area of the outer corona. The signal of these are was 
filtered using the a Trous filtering algorithm and the automated technique of K04 
was used for the detection of intensity oscillations. 1054 pixels found to oscillate in 
intensity. 



out of the 4200 pixels of the sample while K04 found 84 (out of which 66 
were concentrated in a very compact area in the middle of the image) . 

The difference in number of detections before and after the is signif- 
icant. The number of detections before and after the a Trous wavelet 
transformation increased by a factor of 75, 96 and 3.3 for the Moon, 
outer corona and lower corona areas respectively. What is also impor- 
tant is that the increase is not the same for the three areas. While 
the increase in number of detections is similar for the Moon and outer 
corona areas, it is by far smaller for the lower corona. 

Although it is not surprising that the a Trous filtering effects the 
areas with high S/N ratio less than those with very low (or zero) S/N, 
it might appear strange that the filtering causes the areas with very 
low S/N ratio to be detected as oscillating. To examine this difference 
in some detail the wavelet transformation of two pixels, one from the 
Moon area and another from the lower corona are included as figures 
5 and 6. Both points were not found to oscillate before the filtering 
but only after. By examining the time series on Panel (a) of Figures 
5 and 6, it is obvious that they are both very noisy, although in the 
case of Figure 6 there is an underling longer timescale variation, while 
in Figure 5 the signal oscillates around an average value. These dif- 
ferences are also apparent in the panel (b) of the two figures were the 
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Figure 4- A 70x60 area of the lower corona and the Moon limb. The area was 
treated with the same algorithm as Figures 2 and 3 and 276 detections were found. 



wavelet transformation in Figure 5 has a lot of high values in very 
low periodicities (since those are far more effected by non-systematic 
noise), while Figure 6 having higher S/N ratio is less effected by noise 
therefor there are less high values in the wavelet transformation even 
in low periodicities. On high periodicities, although there is an area 
of interest in both figures, there is nothing that satisfies the criteria 
established by K03. 

After filtering with the a Trous wavelet transformation algorithm the 
same two points were analysed using wavelets. Figure 7 contains the 
time series produced and and wavelet transformation that corresponds 
to point of the Moon area that we analysed in Figure 5. All the jittering 
in the time series has disappeared and only some small picks and a 
small long-term variation have remained. The wavelet transformation 
corresponds well to what appears on the time series, producing very 
low values on the very high frequencies (as the Gaussian noise effects 
the high frequencies more), some short-lived high values on the high 
frequencies (that correspond to the high picks of the time series) and 
a number of detections in lower frequencies (that correspond to the 
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Figure 5. Wavelet analysis of pixel x=344, y=281 of the Moon area. Panel (a) 
contains the time series while panel (b) the wavelet transformation. This is the 
un-filtered time series and there are no detections of oscillations that satisfied the 
criteria established by K03. 



long-term variation). As expected, the d Trous wavelet transformation 
filtering was very effective on removing the Gaussian noise (which is 
why there are no oscillation in very low periodicities) and the detections 
on the higher periodicities should be attributed to another factor. Since 
by definition the area of the images covered by the Moon has not 
direct light from the lower corona, another source of light should be 
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Figure 6. Wavelet analysis of pixel x=374, y=221 of the lower coronal area. Simi- 
larly to Figure 5, this is the un- filtered time series and there are no detections of 
oscillations. 



considered. As it is known that in total solar eclipses the sky is not 
completely dark (the sky is much brighter that during night even to 
the naked eye), the long-term variations in the time series and the 
resulting detections should be attributed to the scattered light and 
atmospheric affects that produces variations in brightness. Although 
those existed previously in the un-filtered data, there were small and 
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Figure 7. Wavelet analysis of pixel x=344, y=281 of the Moon area after the time se- 
ries was filtered through the a Trous wavelet transformation. A number of oscillation 
can now be found. 



effect by Gaussian noise, therefor they did not produce enough power 
to become valid detections. 

Figure 8 contains the time series and wavelet transform of the point 
of lower corona analysed in Figure 6. The effects of the a Trous wavelet 
filtering here are different to those on the previous time series. Although 
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Figure 8. Wavelet analysis of pixel x=374, y=221 of the lower coronal area after 
the time series was filtered through the d Trous transformation. As with Figure 7, 
a number of oscillations can be found. 

most of the jittering has disappeared up to the 130 th sec, the variation 
that existed on the un-filtered data after the 130 th sec still largely 
remains. The same applies to the transformation of panel 9b) of Figure 
8: All high values in very low periodicities have disappeared up to 130 th 
sec, but a significant number remains after then. On higher periodicities 
there is clearly an amplification of existing high values mainly on the 
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region after the 130 sec. The increase in power on the low periodicities 
coincides with a similar increase in the power of high periodicities. 



4. Significance test using Monte Carlo analysis 

The statistical analysis on SECIS 2001 as performed by K03 & K04 
indicates that there are "false" (i.e. caused by noise or atmospheric 
effects) oscillations that satisfy the criteria establish by K03. In par- 
ticular atmospheric seeing is known to produce differential distortion 
effects that produce "false" detections in the range of 5 to 19 sec 
periods. So far the only satisfactory way found to treat those false 
detections is statistical. The possibility for a detection to be due to 
noise or atmospheric effects was calculated by scanning for oscillations 
large parts of the data were no "real" (i.e. solar coronal) oscillations 
can be expected. The number of oscillations found in those areas was 
then used to establish the possibility of detection to be false in the 
areas of the lower solar corona that are in the proximity of the "test" 
areas. The cases were the concentration of detections is much higher 
than expected, a detection of a corona wave is reported (as in K04). 

In a bid to establish a quantitative criterion to determine which 
detections of oscillations are due to noise or atmospheric effects and 
overcome the limitations of the statistical methods used, the Monte 
Carlo analysis (or randomisation) was investigated. This method has 
successfully used before in the analysis of solar physics data by O'Shea 
et al., Banerjee (2001, 2001) who applied it to time series analysed 
by wavelet transforms. It was chosen because it provides a test of 
noise that is distribution free (or else non-parametric), i.e. it does 
not depend on any given noise model (e.g. Gaussian noise, Poisson, 
etc). Here we follow the Fisher method as described by Nemec & Ne- 
mec (1985) and performed 1000 permutations per pixel of the aligned 
three-dimensional array. In order to evaluate the performance of this 
randomisation method, we used the same test areas of the Moon and the 
outer corona as in previous section. For each individual pixel of these 
areas the maximum power of the wavelet transformation was recorded 
and then compared to the maximum value of the 1000 shuffled time 
series produced from the original. The percentage p of those shuffled 
time series that had maximum power of their wavelet transformation 
larger than that of the original time series was then recorded. When 
the original time series is random noise of any given type (as this test 
is distribution free) we expect that 50% of the shuffled time series will 
have wavelet transformations with higher maximum power than the 
original data. For the purposes of this analysis we will consider any 
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Table I. Percentage of pixels in each of the test areas before and after the 
noise filtering that were found to have p > 1%. 



Percentage of p > 1% (%) Before noise filtering After noise filtering 



Moon 
Outer Corona 
Lower Corona 



60 
65 
5 



15 
10 
1 



value p smaller than 1% (ie less than 1% of the shuffled time series has 
wavelet transformation with power more than that of the original time 
series) is indicating that the original time series has a strong signal 
comparing to the noise level. 

Table I contains the percentage of pixels in each of the test areas 
that were found to have p > 1%. A number of important results become 
apparent: First of all the areas of the Moon and outer corona before 
noise-filtering have more that the half of their pixels having a relatively 
low S/N ratio, while the percentage of the lower corona pixels (again 
before filtering for noise) that have a similarly low S/N ratio is much 
lower (5 %). This is significant as it confirms the choice of p > 1% as a 
criterion to distinguish between the pixels that have signal dominated 
by the solar corona and those that do not. Second the application of 
the d Trous noise filtering makes reduces the "randomness" of the data 
set significantly (by a factor of ~ 5) making it difficult to use the 
randomisation test reliably. A third useful result shown on Table I is 
that the percentage of pixels that have p > 1% is approximately the 
same in the Moon area (where we know that the signal of all pixels 
is due to not direct observation of the Sun) and the outer corona 
(where we know the signal is partially directly from the solar corona 
and partially scattered light). This is important as it indicates that the 
scattered light on this image area is a significant portion of the signal 
and it will be very difficult to distinguish which detections are "real" 
in this area by applying the existing criteria. 

A Trous wavelet transformation's inability to reduce the number of 
"false" detections can be explained if we consider the different contri- 
butions to the signal. Any pixel value of this data set is a combination 
of the detection of scattered light from Earth's atmosphere, Gaussian 
noise and (for those areas that are not covered by Moon) the detection 
of light directly from the solar corona. By removing the Gaussian noise 
any weak oscillations due to Earth's atmosphere could be detected more 
clearly. Those oscillations are known to be introduced by various opti- 
cal effects produced by Earth's atmosphere: variation in transmission 
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through the atmosphere, differential distortions caused by winds at high 
altitude, etc. In contrast the area lead by the lower solar corona had 
already enough signal and a high S/N ratio, therefor the reduction of 
the Gaussian noise levels did not contribute to a major increase in the 
number of detected oscillations. Also because the signal coming from 
the solar corona in this particular area was strong before the a Trous 
noise filtering, the relatively weaker atmospheric effects did not increase 
dramatically after the subtraction of the Gaussian noise. As a result 
the number of detected oscillations did not increase as dramatically in 
that area as in the other two. 

The most significant breakthrough of the efforts to establish a quan- 
titative criterion to determine which detections of oscillations are "real" 
and which are not, came when p was calculated for the pixels that were 
found to oscillate by K04 (i.e. the 5 pixels of the Moon area, 11 pixels 
of the outer corona area and 66 pixels of the lower corona). All pixels 
from lower corona were found to have p < 0.1% (i.e. none of the 1000 
shuffled time series was found to have a wavelet transformation with 
higher power than that of the original time series), while 4 out of 5 of 
the detections from the Moon area and 10 out of 11 detections from the 
outer corona, were found with p > 0.1%. Therefor, a criterion can be 
establish that will use the randomisation test described here to reject 
those pixels with p > 0.1%. 



5. Discussion 

Two well known signal processing techniques, The a Trous wavelet 
transformation and Monte Carlo analysis, were applied to SECIS 2001 
data. Those two methods were evaluated by using two "test" areas 
(areas were the signal from the solar corona was expected to be small 
or zero) and a "useful" area were detections of corona oscillations were 
expected. By comparing the results from the three areas, an accurate 
evaluation of the numerical techniques described above was made. 

The a Trous algorithm produced mixed results. Although the reduc- 
tion to Gaussian noise level was very significant, the ability to detect 
corona waves was actually reduced. This is because of the effect of 
Earth's atmosphere in the data set. Intensity oscillations caused by 
the atmosphere were weaker in signal than those caused by the solar 
corona, therefor, when the S/N ration was lower those oscillations did 
not obtain the significance levels needed to become detections. After 
the noise was reduced the significance levels of the intensity oscillations 
caused by the atmosphere was increased enough to produce a large 
number of false detections. It is also worth noticing that the areas where 
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the signal from the corona was stronger had much less false detections 
than those areas with weak or no signal, indicating that the signal from 
the lower corona is significantly stronger than the atmospheric effects. 

By the use of the two "test" areas of the data set it has become 
apparent that a reliable, objective, numerical method is needed in order 
to distinguish those detections caused by plasma from the solar corona 
to those introduced by atmospheric effects. The Monte Carlo analysis 
(otherwise refereed here as randomisation test) was investigated as a 
means to make this distinction. All detections reported as "real" by K04 
were found to have p < 0.1% while almost all (14 out of 16) of those 
reported as "false" were found to be in the range of 100% > p > 0.1%. 
Therefor the value of p is proposed as a criterion for rejecting future 
detections from the lower corona region. 
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